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The  estimated  residual  error  (or  its  bound)  for  numerical  quadratures  is  usually  expressed  in  terms  of  a  derivative  of 
some  order  of  the  integrand  or  some  residual  factor  of  the  integrand  after  factoring  out  a  countable  number  of  zeroes  and 
singularities  that  occur  along  the  integration  path.  The  order  of  the  derivative  is  a  function  of  the  number  of  sample 
points  for  evaluating  the  integrand.  Regrettably,  the  magnitudes  of  these  higher  order  derivatives  are  difficult  enough  to 
estimate  for  even  analytic  sound  velocity  profiles.  In  practice,  observed  sound  velocity  profiles,  which  are  usually  given  in 
tabular  form  and  include  measurement  errors,  exacerbate  our  inability  to  assess  the  magnitudes  of  these  higher  order 
derivatives.  An  estimate  of  the  residual  error  expressed  in  terms  of  a  first  derivative  would  be  far  more  practical  for  both 
analytic  and  observed  sound  velocity  profiles.  f 
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CONVERGENCE  RATE  OF  CODES 


FOR  NUMERICAL  QUADRATURE  TECHNIQUES 
FOR  CLASSICAL  RAY  TRACING 

Edward  R.  Floyd 
Arctic  Submarine  Laboratory 


Range  Integral : 


C  ( 2 ) 

CcS  -  c*(z)]i/z  dz 

pv 

=  I  W(z)  f(z)  dz 

Vi  and  f  are  convenient  functions  of  z.  [W  dz]  becomes 
a  Stieltjes  measure  and  f(z),  the  Stieltjes  integrand. 

Trick:  Choose  W(z)  to  incorporate  any  poor  behavior 
(singularities  or  zeros)  into  W  leaving  f(z)  well 
behaved. 


Thus  the  Gaussian  quadrature 


n  is  the  order  of  the  approximation; 

m  is  determined  by  W; 

Wj  n  is  a  Gaussian  weighting  coefficient 
determined  by  W; 

n  is  a  Gaussian  sample  point  determined 

by  VL-. 


The  Gaussian  error  is  given  by 
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where  z  <  K  < 
o  —  — 


For  underwater  acoustics 


Factors  cut  branch  point  singularity  at  vertex 


Therefore 
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This  expression  has  several  deficiencies; 


*  Usually  impractical  to  evaluate  f 


(2n) 


f^1*)  for  i  <  2n  may  be  discontinuous, 


C(z)  is  often  known  only  in  tabular  form 
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Analytic  (hypothetical)  C(z): 


•  Rapid  convergence  with  increasing  n. 

•  Only  limited  by  round-off  error. 

•  Thumb  rule  for  maximum  effective  n: 

nmax  ~  ("number  of  signifcant  figures" )/2 


E.  R.  Floyd,  J.  Acoust.  Soc.  Am.  49,  1580-1590 
(1971) 
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Observed  (nonanalytic)  C{z): 


•  C(z)  tabled  at  discrete  depths. 

’  Cubic  spline  fit  — ►  d^C/dz^  discontinuous. 
Nevertheless  still  rapidly  convergent  over  an 
assembly  of  116  observed  C(z)'s.  MOE  —  t\ 
reduced  by  a  factor  of  0.763  for  each  additional 
sample  point  per  Riemann  sheet. 

*  Thumb  rule  for  maximum  effective  n: 

nmax  :L "number  of  significant  figures". 


*  E.  R.  Floyd  and  J.  D.  Pugh,  J.  Acoust.  Soc.  Am. 
61  ,  682-687  (1977). 
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Present  Status: 


Even  for  analytic  functions,  quadrature 
errors  are  difficult  to  estimate  because 
higher  order  derivitives  are  generally 
comp! icated. 

While  we  have  convergence  rates  for  an 
assembly  of  observed  profiles,  the 
discontinuity  of  higher  order  derivatives 
confounds  our  estimate  of  convergence  for  a 
particular  observed  profile. 

Desire  an  estimate  expressed  in  terms  of  the 
first  derivative  of  the  Stieltjes  integrand. 


